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I.  INTRODUCTION 

The  energy-transfer  mechanisms  for  self-deactivation  of  vibrationally 
excited  hydrogen  fluoride  molecules  are  quite  complex.  Possible  mechanisms 
include  energy-transfer  processes  in  which  (1)  vibrational  energy  is  exchanged 
between  colliding  HF  molecules  (v  -*  v)  through  intermolecular  energy  ex- 
change; (2)  vibrational  energy  is  converted  into  rotational  and  translational 
energy  (v-»  R,  T)  through  intermolecular  or  intramolecular  energy  exchange, 
or  both;  (3)  rotational  energy  is  converted  into  vibrational  energy  (R  -*  v) 
through  intramolecular  energy  exchange;  and  (4)  rotational  energy  is  con- 
verted into  rotational  and  translational  energy  (R  -»R,  T)  through  inter- 
molecular and  intramolecular  energy  exchange.  It  is  difficult,  if  not  impos- 
sible, to  determine  experimentally  detailed  rate  coefficients  for  these  four 
possible  types  of  energy- transfer  mechanisms.  A knowledge  of  the 
temperature-dependent  rate  coefficients  for  these  types  of  energy-transfer 
mechanisms  is  important  for  an  understanding  of  HF  chemical  lasers.  Some 
of  the  problems  associated  with  interpretation  of  previous  work  on  the  self- 
deactivation  of  vibrationally  excited  HF  are  briefly  summarized  here. 

A.  'HF(vt  = 1)  + HF(v2  = 0)  SYSTEM 

Vibrational  relaxation  rates  of  HF(vj  = 1)  by  HF(v£  = 0)  molecules  have 

been  measured  by  several  techniques.  * The  deactivation  rate  of  HF(Vj  = 1)  by 

12  3 -1  -1 

HF(v2  = 0)  has  an  average  measured  value  of  1.2  X 10  cm  mole  sec  at 

2 

room  temperature.  In  the  interpretation  of  the  experiments  reported  in  the 


I 

I 


published  literature,  it  is  assumed  by  most  authors  in  reducing  their  data 

i i 

that  the  reaction  HF(vj  = 1)  + HF(v£  = 0)  -*  HF(vj  = 0)  + HFfv^  = 0)  occurs 

through  a v T,  R energy- transfer  process,  which  has  an  energy 

defect  of  3962  cm'*;  this  energy  is  assumed  to  be  given  up  to  the  thermal 

3 

bath.  Airey  and  Fried  postulated  that  these  fast  self- relaxation  rates  of 
HF  are  the  result  of  a strong  attraction  potential  that  exists  between  HF 
molecules,  which  results  in  the  formation  of  HF  dimers,  and  that  they, 

3 

perhaps,  could  be  caused  by  monomer-dimer  collisions.  Airey  and  Fried 
point  out,  however,  that,  at  the  pressures  and  temperature  of  their  experi- 
ment, less  than  0.001%  of  the  HF  exists  as  dimers,  so  that  the  measured 

rate,  in  their  opinion,  is  the  result  of  monomer -monomer  collisions. 

4 

Thompson  points  out  that  Airey  and  Fried’  s point  of  view  is  correct 
if  it  is  assumed  that  the  vibrational  re laxati on. occurs  by  collisions  between 
the  excited  HF  molecule  and  the  HF  dimer,  but  it  is  not  necessarily  correct 
if  the  effective  intermediate  is  a long-lived  complex  formed  by  collision  of 
an  excited  HF  molecule  with  another  HF  molecule.  Under  the  latter  cir- 
cumstances, the  vibrational  relaxation  rates  are  critically  dependent  on  the 
rate  of  formation  and  dissociation  of  the  dimer.  The  long-lived  complex 
corresponds  to  two  HF  molecules  held  together  by  an  attractive  well.  The 
complex  can  be  stabilized  temporarily  by  transfer  of  collision  energy  into 
several  modes  of  the  complex  such  as  rotation.  Transfer  of  energy  from 
these  modes  back  to  translation  results  in  an  equilibrium  condition  between 
the  monomers  and  dimers  of  HF. 
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Shin  used  a vibratidn-to- rotation,  energy-transfer,  one-dimensional 
analytical  model,  which  included  the  effects  of  the  dipole-dipole  and  hydrogen- 
bond  iftteractions  to  explain  the  efficient  vibrational  relaxation  of  HF(v^  = 1) 
by  HF(v£  = 0)  molecules  at  low  temperatures.  Shin  speculated  that,  because 
of  the  strong  dipole -dipole  attractive  forces,  the  HF  molecules  attract  each 
other  at  large  separations  even  when  each  HF  molecule  is  rotating  rapidly. 

The  addition  of  the  dipole-dipole  interaction  term  to  the  short-range  inter- 
action term  enhances  the  probability  of  energy  transfer.  In  HF-HF  collisions. 
Shin’s  model  predicts  that  the  v -*  T energy-transfer  process  is  less  efficient 
than  the  v-»  R energy-transfer  process.  In  the  Shin  analytical  model,  the 
dipole-dipole  and  hydrogen-bond  interactions  are  considered  as  separate, 
important  interactions  in  the  relaxation  of  vibrationally  excited  HF. 

One  objective  of  this  research  was  to  investigate  the  requirement  that 
HF  dimers  be  formed  at  room  temperature  in  order  to  explain  the  fast 
self-relaxation  rates  measured  by  various  experimental  techniques.  * A 
second  objective  was  to  calculate  the  temperature-dependent  rate  coefficient 
for  the  vibrational  relaxation  of  HF(v^  = 1)  by  HFfv^  = 0)  and  compare  the 
result  with  available  experimental  data.  A potential  energy  surface  for 
the  collisions  between  HF(vj,  J^)  + HF(v£,  J^)  molecules  was  constructed, 
which  consists  of  a London- Eyring-Polanyi-Sato  (LEPS)  potential -energy 
function  for  the  short-  range  interactions  and  a partial-point  charge,  dipole - 
dipole  function  for  the  long-range  interactions  between  the  four  particles. 
Three-dimensional  trajectory  calculations  were  made  of  the  collision 
dynamics  of  the  HF(Vj  = 1,  Jj)  + HFfv^  = 0,  3^)  systems  on  this  energy 
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surface.  Temperature -dependent  rate  coefficients  are  reported  for  the 
v-*  R,  T energy- transfer  processes  that  occur  in  collisions  involving 
HF(Vj  = 1,  Jj)  + HF(v2  = 0,  J^)  systems. 

B.  HF(Vl  > 1)  + HF(v2  = 0)  SYSTEMS 

Knowledge  of  the  temperature-dependent  rate  coefficient  for  vibrational 
relaxation  of  HF(Vj  = 1)  by  HF(v2  = 0)  is  not  sufficient  to  model  the  HF  chem- 
ical laser.  It  is  necessary  to  know  also  temperature -dependent  rate  coeffi- 
cients for  vibrational  relaxation  of  HF  species  from  the  higher  vibrational 
levels  (Vj  = 2 through  Vj  = 6).  Rate  coefficients  that  include  contributions 
from  both  v -*  v and  v -*  R,  T energy- transfer  processes  for  the  deactivation 
from  the  higher  vibrational  levels  of  HF  have  been  reported  by  several 
authors.  ~ Some  of  these  authors  ’ estimated  the  contribution  to  the  rate 
coefficients  for  the  energy-transfer  processes  by  means  of  a harmonic 
oscillator  model  k . = vk4  n to  predict  the  contribution  to  the  total  rate 
coefficient  from  v-*-  R,  T energy  transfer.  In  the  harmonic  oscillator  model, 
the  v -•  R,  T energy- transfer  processes  do  not  involve  multiquantum  v -*  R,  T 
transitions.  The  rate  coefficients  for  the  endothermic  v-»  v energy- transfer 

i- 

processes  HF(Vj  > 1)  + HF^  = 0)  -*  HF(vj  = Vj  - 1)  + HF(v'  = 1)  can  be  cal- 
culated if  the  rate  coefficients  are  available  for  the  v -»  R,  T energy-transfer 

8 

processes.  Flowtube  measurements  reported  by  Kwok  and  Wilkins  indicate 
2 

a v dependence  for  the  rate  coefficients  that  represent  the  sum  of  all 
v -»  R,  T energy  processes  that  occur  from  a given  v level  for  the  DF(v^)  + 

HF(v2  = 0)  system.  Since  the  dipole-dipole  interactions  between  the  colliding 
HF  + HF,  DF  + DF,  and  the  HF  + DF  systems  are  quite  similar,  it  is 

1 
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probable  that  the  v -*  R,  T total  rate  coefficients  for  DF(Vj)  + DF(v2  = 0),  and 
HF(vj)  + HF(v2  = 0)  exhibit  a v dependence  similar  to  that  for  DF(v^)  + 

HF(v2  = 0). 

It  has  not  been  possible  to  separate  the  contributions  from  the  v-.v  and 
v -»  R,  T energy-transfer  processes  obtained  in  experimental  measurements. 

aje  * 

In  an  HF  + HF  system,  where  HF  represents  vibrationally  excited  HF 
* 

species,  the  HF  vibrational  levels  may  relax  at  the  same  time.  If  the  v-*  v 

rates  are  faster  than  the  v-»  R,  T rates,  the  v -*  v rates  would  couple  together 

all  the  vibrational  levels  so  that  the  rate  coefficients  for  the  v -*  R,  T 

energy-transfer  processes  could  not  be  separated  from  the  rate  coefficients 

for  the  v -•  v energy-transfer  processes.  Only  the  kj  ^ rate  coefficient  for 

v-»  R energy  transfer  can  be  obtained  without  this  v _*  v coupling  difficulty. 

There  is  a definite  need  for  both  the  temperature  and  v dependences  of 

the  rate  coefficients  for  both  v -»  v and  v-*R,  T energy-transfer  processes. 

Experimental  v -*  v rate  coefficients  can  be  deduced  from  measured  (v  _ v + 

v -»  R,  T)  rate  coefficients  only  when  information  is  available  on  the  v-»R,  T 

2 

rate  behavior  of  the  upper  v levels.  A v or  greater  dependence  in  the  v -*  R 
rate  coefficients  would  result  in  a further  decrease  in  the  v-»  v rate  coefficients 
and  would  probably  eliminate  the  problem  of  predicting  from  experiment 
v -*  v rate  coefficients  that  exceed  the  gaskinetic  collision  rate.  The  com- 
puter modeling  of  flowtube  results  of  HF(v^  > 1)  + HF(v2  = 0)  systems  pro- 
vides indirect  evidence  that-  such  large  v-*v  rate  coefficients  for  the  upper 
vibrational  levels  are  not  correct.  When  such  large  v-»  v rate  coefficients 
are  used  to  model  the  flowtube  results,  it  is  not  possible  to  obtain  any 


-13- 


agreement  with  the  measured  decay  slopes.  Additional  experimental 
evidence  that  such  large  v-*  v rate  coefficients  are  not  correct  involves  the 
measurement  of  the  time- resolved  spectral  behavior  for  the  v -*  v - 1 
transitions  at  room  temperature.  Too  much  laser  energy  output  is  pre- 
dicted from  the  upper  vibrational  levels  when  such  large  rate  coefficients 
for  the  v-*  v rates  are  used  in  pulsed  or  continuous-wave  chemical-laser 
computer  codes.  The  effect  of  fast  v -*  v rates  is  to  pump-up  the  higher 
vibrational  levels  more  than  would  be  predicted  on  the  basis  of  the  initial 
product  distributions  measured  for  the  pumping  reactions  used  to  generate 
the  lasing  species.  It  would  seem,  therefore,  that  the  v -*  R,  T rate 
coefficients  must  have  a greater  than  the  v dependence  predicted  from  the 
harmonic  oscillator  model.  The  . v rate  coefficients  for  the  exothermic 
processes  HF(vj  > 1)  + HF(v£  = 1)  -«HF(Vj  = Vj  + 1)  + HFfv^  = 0)  would  not 
be  larger  than  the  gaskinetic  collision  rates,  as  predicted  from  recent 
f xperimental  measurements,  if  (1)  the  v -»  R,  T rate  coefficients  for  the 

endothermic  processes  HF{v^  > 1)  + HF(v£  = 2)  -*  HF(vj  - 1)  + HF(v£  = 1) 

2 

have  a v or  greater  v dependence  and  (2)  if  the  energy  mismatch  AE  used 
to  calculate  the  rate  in  the  exothermic  direction  includes  the  effect  of 
rotational  motion  of  the  colliding  HF  molecules. 

A one -dimensional  theory  of  vibrational-to-vibrational  energy  transfer 
was  developed  by  Shin.  In  this  theory,  it  is  assumed  that  energy  transfer 
occurs  through  the  formation  of  nonrigid  dimers  at  low  temperatures 
(T  £ 300  K)  and  through  rotational  motion  of  the  colliding  molecules  at  high 
temperatures  (T  > 300  K).  In  the  nonrigid  dimer  model,  the  energy  mismatch 
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is  transferred  to  the  restricted  translation  motion  of  die  colliding  molecules, 
and,  in  the  rotational  model,  the  energy  mismatch  is  transferred  to  die 


rotational  motion  of  the  colliding  molecules.  The  rate  coefficients  predicted 

10  » 
by  Shin  for  the  deactivation  processes  HF(vj)  + HF(v£  = 0)  -*  HF(vj  = 

t 

Vj  - 1)  + HF(v£  = 1)  at  room  temperature  increase  with  v^,  but,  after  reaching 
a maximum  value  at  v^  = 4,  they  decrease  sharply  with  increasing  v^.  Shin 
and  Kim**  extended  Shin's  one- dimensional  model***  to  study  the  temperature 
dependence  of  the  energy- transfer  rates  in  die  temperature  range  200  to 
2000  K.  They  found  that;  at  200  K,  the  rate  coefficients  for  those  v-*  v pro- 
cesses decrease  with  increasing  v because  of  the  dimer  mechanism;  at  200  K, 
a maximum  value  in  these  rate  coefficients  occurs  at  v^  = 3;  and,  above  300  K 
the  rate  coefficients  increase  with  increasing  Vj  as  a result  of  the  rotational 
contribution  of  the  rotational  motion  to  the  energy-transfer  processes.  Shin 
and  Kim's  model,  when  applied  to  two-quantum  v-*v  processes  HF(vj)  + 

i i 

HF(v£  = 0)  -» HF(vj  - 2)  + HFfvg  = 2),  predicts  very  small  values  of  the  rate 
coefficients.  Their  theory  does  not  indicate  the  mechanism  by  which  vibra- 
tional energy  is  converted  into  rotational  energy  in  HF-HF  collisions,  nor 
does  it  predict,  for  HF.  HF  collisions,  multiquantum  vibrational- rotational 
energy  exchange  processes. 

One  objective  of  this  research  was  to  investigate  the  mechanism  by 
which  the  vibrational  energy  of  the  vibrationally  excited  HF  molecule  is 
converted  into  rotational, energy.  Another  objective  was  to  calculate  the 
temperature  and  v dependences  of  the  rate  coefficients  for  both  the  v -*  v 
and  v-*  R,  T energy-transfer  processes  in  HF(v^  > 1)  + HF(v^  = 0)  collisions. 
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C.  ROTATIONAL  RELAXATION  OF  HFfVj,  Jj)  BY  HF(v2  = 0,  J.J 

12  13 

Green  and  Hancock  ’ measured  HF(v  = 1)  self- relaxation  rates  in 
pure  HF  and  in  HF  highly  diluted  with  Ar  by  means  of  the  laser-excited 
vibrational  fluorescence  technique.  They  found  that  the  HF(v  = 1)  self- 
relaxation rate  observed  in  pure  HF  was  30%  slower  than  the  similar  rate 
with  the  HF  diluted  in  Ar.  They  gave  one  explanation  of  this  behavior  as  being 
the  result  of  incomplete  rotational  thermalization  of  high  rotational  J levels 
during  vibrational  relaxation  in  pure  HF.  With  it  assumed  that  a quantum  of 
vibrational  energy  of  the  HF(v  = 1)  molecule  goes  preferentially  into  rotational 
energy  of  the  same  molecule  upon  an  HF(vj  = 1)  + HF(v£  = 0)  collision,  such 

an  energy  mechanism  would  give  a large  AJ  change  with  a very  small  energy 

i r 

defect.  The  rotationally  excited  HF(Vj  = 0,  Jj  = 14)  molecule  produced  can 
be  expected  to  have  a lifetime  comparable  to  the  vibrationally  excited 
HF(vj  = 1)  molecule  vibrational-relaxation  time.  If,  for  example,  the 

i i 

HF(vj  =0,  Jj  = 14)  molecules  produced  by  the  v -*  R mechanism  are  not 

removed  immediately,  it  is  then  possible  for  them  to  transfer  back  to 

HF(Vj  = 1,  J = 2)  by  an  R -*  v process  (tile  result  of  small  energy  defects, 

i.e. , less  than  200  cm”*)  and  consequently  decrease  the  apparent  measured 

HF(vj  = 1)  self- relaxation  rate  by  the  v-*  R,  T energy-transfer  processes. 

14 

This  mechanism  is  not  as  effective  when  Ar  is  used,  since  Ar  efficiently 

rotationally  quenches  HF  and  DF. 

14 

Hancock  and  Saunders  studied  the  self- relaxation  of  DF(v  = 1)  and 
found  that  their  room- temperature  results  clearly  indicate  a 50% 


slower  DF(v  = 1)  self-quenching  rate  in  pure  DF  than  when  measured  in 

DF-Ar  mixtures.  This  decrease  in  the  relaxation  time  of  DF(v  = 1)  in  pure 

DF  can  be  explained  if  it  is  assumed  that  the  vibrational  and  rotational 

relaxation  occur  on  similar  time  scales  in  undiluted  DF.  Hancock  and 
13 

Green  speculated  that  Ar  efficiently  quenches  rotationally  excited  HF  and 

DF.  This  means  of  rotational  deactivation  of  HF(v  = °'  Jhigh> or  DF(v  * °- 

Jhigfa)  with  Ar  is  not  applicable  when  studying  vibrational  relaxation  of 

HF(v  = 1)  or  DF(v  = 1)  in  pure  HF  or  DF,  respectively.  The  complicated 

situation  that  the  vibrational  relaxation  and  the  rotational  relaxation  bccur 

on  similar  time  scales  must  be  dealt  with. 

14 

Hancock  and  Saunders  reported  that  their  measured  HF(vj  - 1)  - 

DF(v2  = 0)  quenching  rates  indicated  no  discrepancy  between  the  diluent 

versus  no  diluent  cases.  The  argument  used  for  pure  HF  or  DF  can  be  used 

also  to  explain  the  lack  of  difference  observed  in  their  HF(vj  = 1)  + DF(v2  = 0) 

quenching  rates.  Most  of  the  vibrational  energy  in  HF(v^  = 1)  + DF(v2  = 0) 

15  1 

collisions  is  transferred  through  v -*  v exchange,  ' with  an  energy  defect 
of  1055  cm-*.  The  v-»  R energy- transfer  processes  are  still  important, 
but  the  energy  that  must  be  absorbed  [1055  cm-*  for  HF(vj  = 1)  + DF(v  = 0) 
compared  to  3962  cm-*  for  HF(v^  = 1)  + HF(v2  = 0)  and  2907  cm-*  for 
HF(vj  = 1)  + DF(v2  = 0)]  is  much  smaller.  Lower  rotational  states  of 
DF(vj  = 0)  would  be  populated  and  the  rotational  energy  transfer  would 
proceed  faster  because  of  the  close  spacings  of  the  low  rotational  energy 
levels.  The  rate  coefficient  for  the  deactivation  of  DF(v^  = 1)  by  DF(v2  = 0) 
is  approximately  less  than  half  the  rate  coefficient  for  deactivation  of 
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HF(v  = 1)  by  HF(v  = 0),  ^ which  is  additional  evidence  that  nearly  all 

of  the  initial  vibrational  energy  of  HF(v  = 1)  is  converted  into  rotational 

energy  rather  than  translational  energy  of  the  colliding  pair.  This  result 

concurs  with  earlier  results  for  HC1  and  DC1  presented  by  Chen  and 
17 

Moore. 

The  rotational  relaxation  of  highly  rotational  excited  HF  molecules  has 

18 

only  been  investigated  qualitatively.  Polanyi  and  Woodall  and  Ding  and 
19 

Polanyi  showed  that,  for  hydrogen  halides,  the  probability  of  rotational 
deactivation  decreases  with  increasing  rotational  quantum  number  J.  The 
vibrational  energy  of  vibrationally  excited  HF  molecules  can  be  converted 
almost  entirely  into  rotational  energy  by  a mechanism  that  populates  the 
high  rotational  quantum  states  of  the  same  vibrationally  excited  HF  molecules 
(large  AJ  change)  with  small  energy  defects.  It  is  expected,  therefore,  that 
the  rate  coefficients  for  deactivation  from  these  high  J states  by  R — T 
processes  decrease  with  increasing  J because  an  increase  in  J is  associated 
with  an  increase  in  the  energy  spacings  between  rotational  levels.  In  addition 
to  R -*  T energy  transfer,  it  is  expected  that  collisions  with  resonance 
transfer  of  rotational  energy  play  an  important  role  in  rotational  deactivation. 
This  type  of  collision  is  R -*  R energy  transfer.  The  rotational  energy  given 
up  by  one  HF  molecule  is  transferred  to  its  partner  in  collision.  The  R -*  R,T 
energy-transfer  processes  provide  the  drainage  necessary  for  the  rotational 
nonequilibrium  states  to  relax  to  a thermal  distribution. 


-18 


.WM..,.  .,L.,,  — « - » wwmxm mm 


20 

Deutsch  observed  laser  action  on  pure  rotational  transitions  of  HF 

formed  by  the  chemical  reaction  F + H£.  He  observed  rotational  transitions 

in  HF  from  energy  levels  corresponding  to  very  high  J values,  indicating  a 

21 

rotational  nonequilibrium  distribution  of  HF(v,  J)  states.  Polanyi  points 

out  that  this  rotational  nonequilibrium  can  exist  when  the  spacing  between 

the  rotational  levels  becomes  large  enough  such  that  a few  collisions  are  not 

sufficient  to  induce  rotational  relaxation  by  AJ  = 1 rotational  jumps,  as  is 

20 

the  case  when  the  rotational  spacing  is  small.  Deutsch  reports  rotational 

lasing,  for  example,  from  v = 0,  J = 28,  from  v = 2,  J = 20,  and  from  v = 3, 

22 

J = 28.  Theoretical  trajectory  calculations  by  Wilkins  for  the  F + H^  reaction 

ii  ii 

predict  that  the  HF(v  , J ) states  are  populated  up  to  J =17  for  v = 1 and 

' i 23 

up  to  J = 12  for  v = 2 or  3.  Polanyi  and  Tardy,  in  their  study  of  the 
F + Hg  reaction  by  means  of  an  infrared  chemiluminescence  technique, 

i i 

observed  HF  emission  from  v =1,2,  and  3,  with  J = 12,  12,  and  4, 

’ ' ' max 

i 

respectively.  No  HF  emission  was  observed  from  v = 4.  From  thermo- 
dynamic considerations,  the  reaction  energy  of  the  F + system  is  sufficient 

i i 

to  populate  v =1,2,  and  3 up  to  J = 20,  14,  and  5,  respectively.  The 

rnstx 

rotational  states  responsible  for  the  laser  action  observed  by  Deutsch  could 
not  have  been  formed  from  the  F + H^  reaction.  Another  objective  of  this 
study  was  to  attempt  to  determine  whether  or  not  the  high  J states  are  formed 
by  energy-transfer  processes  involving  vibrational  relaxation  from  higher 
(v,  J)  levels  where  multiquantum  v-*  R processes  occur.  If  multiquantum 
v -*  R energy-transfer-processes  occur,  they  would  provide  a plausible 
explanation  for  the  high  J states  observed  by  Deutsch. 
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Still  another  objective  of  this  report  is  to  give  the  results  of 

calculations  of  temperature- dependent  rate  coefficients  for  the  rotational 

relaxation  of  HF  species  from  both  low  and  high  J states.  In  a subsequent 
24 

report,  these  rate  coefficients  calculated  for  tile  v-*  v and  v -»  R,  T and 
R-*R,  T energy- transfer  processes  will  be  incorporated  into  a model  for 
rotational  nonequilibrium  to  see  if  the  constructed  model  can  explain  some 
of  the  HF  laser  kinetics.  In  Section  U,  the  semiempirical  treatment  of  the 
potential  energy  surface  for  HF-HF  collisions  is  described.  In  Section  III, 
tile  classical  trajectory  calculations  are  discussed.  Results  and  discussions 
of  this  study  are  given  in  Section  IV.  A summary  of  the  results  is  given 
in  Section  V. 


j 
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H.  POTENTIAL  ENERGY  SURFACE 


Before  modeling  die  collision  dynamics,  a reasonable  potential  energy 

surface  must  be  constructed.  The  potential  employed  here  for  HF-HF 

interactions  is  constructed  from  two  functions,  a London- Ey  ring- Polanyi- 

Sato  (LEPS)  potential  energy  function  for  the  short- range  interactions  and  a 

point-charge,  dipole-dipole  potential  energy  function  for  the  long-range 

interactions  between  the  four  atoms.  This  surface  is  capable  of  supporting 

a quasi-bound  state,  i.  e. , the  HF  dimer.  It  has  been  proposed  by  some 
4 5 27-29 

authors  ’ ’ that  the  formation  of  quasi-bound  complexes  in  HF-HF 

collisions  is  responsible  for  the  efficient  v -*  v energy  transfer  observed  in 
several  experiments.  It  is  hoped  that  the  surface  used  in  this  analysis  will 
be  able  to  uncover  the  role  of  dimer  formation  in  HF-HF  collisions 
involving  v-*v,  v-*R,  and  R -*  R energy-transfer  processes. 

The  total  potential  energy  E for  a given  configuration  of  the  four  atoms 
is  defined  by 


where  Ej_,I2ps  *s  potential  energy  function,  and  E^  d is  the  point- 

charge,  dipole-dipole  potential  energy  function.  The  T-EPS  potential  function 
for  the  interaction  between  the  four  atoms  is  defined  by 


21 
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ELEPS  = Q12  + Q34  + Q13  + °14  + Q23  + Q24 


'll  [(J12 


+ J34  " J23  " J24* 


+ (J23  + J24  - J13  - J24)  + (Ji3  + J24  - J12  “ J34 


)4« 


(2) 


where  Q,.  and  J. . are  the  Coulomb  and  exchange  integrals,  respectively. 

J y 

The  Q„s  and  J.jS  are  written  in  terms  of  bonding  and  anti  bonding  energies. 
The  Coulomb  terms  are 


Q..(R. .)  =i  e!*.  (R..)  + (\  j E?.(R..) 

wij  ij  2 ij  ij  \i  + K.^  ijv  ijJ 


(3) 


and  the  exchange  integral  terms  are 


J..(R..)=*  e!*.(R..)  - (1  I-g^) Ea.(R. .) 
y y 2 ijv  ij'  y + k.  .y  i3  iy 


(4) 


where  R..  is  the  internuclear  distance  between  atom  i and  atom  j,  K..  is  the 
ij  ij 

b a 

square  of  the  overlap  integral,  & j is  the  bonding  energy,  and  is  the 
antibonding  energy.  The  bonding  energies  or  the  single-state  energies  are 
represented  by  the  Morse  function 


e!\  = D. . 

lJ  y 


• 2a.  .(R..  — R. .)  -a,  .(R. . — RM 

y y - y . 2e  y y 


(5) 


The  anti  bonding  energies  are  represented  by 


Eij  = Aij 


”2aij(^ij  * ^j)  “aij^ij  ” ^ij^ 


(6) 
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The  diatomic  parameters  required  to  calculate  an  LEPS  surface  for  the 

systems  F + H^,  H + F^»  H + HF,  and  F + HF  were  evaluated  previously  by 
22  26  30  31 

Wilkins.  ' ' ' These  parameters  were  determined  in  previous 

trajectory  studies  in  which  each  surface  was  constructed  for  a specific 

reaction  to  match  either  experimental  activation  energy  for  the  overall 

reaction  or  product  energy  distributions,  or  both.  The  values  of  these 

parameters  are  given  in  Table  I.  In  this  study,  it  was  assumed  that  the 

diatomic  potential  energy  surface  parameters  used  in  the  dynamical  studies 

involving  the  species,  H,  and  F were  transferable  to  the  four-atom  HF-HF 

32 

system.  Previous  trajectory  calculations  by  White  and  Thompson  support 
the  validity  of  the  assumption  that  the  LEPS  equation  diatomic  parameters 
are  transferable. 

The  long-range  potential  E^  ^ is  represented  by  the  Coulomb  electro- 
static function 


1 , 1 
^13  R24 


(7) 


where  ft  is  assigned  a value  of  2.0  X 10”^  esu.  The  function  E^  ^ is  con- 

• o 

structed  as  the  sum  of  four  electrostatic  atom-atom  potential  ± 6 ^/R..  terms; 

y 

each  term  represents  a nonbonded  atom  pair  in  the  HF-HF  system.  The 
2 

term  ± 6 /R..  is  consistent  with  the  classical  model  in  which  the  HF  dipole  is 

considered  to  consist  of  two  opposite  charges,  ft  and  -ft.  The  value  to  be 

33 

used  for  ft  was  taken  from  Be  rend  and  Thommarson,  who  reported  on  a 
planar  classical  trajectory  study  of  vibrational  relaxation  of  HF  and  DF  by 


Table  I.  Parameters  Used  to  Construct  Surface  for  Energy 
Transfer  Between  HF(vj,  Jj)  + HFfv^, 


Parameters 

V k'1 


V k 


D , kcal/mole 
c 


Dipole  Moment 


A..,  kcal/mole 


H-F 

H-H 

F-F 

2.  232 

1.942 

2 

i.  902 

0.  9171 

0.  7419 

1 

L.  4170 

140.  7 


1.  82‘ 


0.  15 


70.  4T 


109.  2 


0.  15 
54.  6^ 


38.02 

0 

0.05 
46.  0C 


aR.  Weiss,  Phys.  Rev.  131,  659  (1963). 

**A.j  = 1/2  De  for  antibonding  states  of  HF  and  H^. 

CR.  L.  Wilkins,  J.  Chem.  Phys.  58,  2326  (1973). 
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HF  and  DF.  These  authors  determined  6 for  the  potential  , by  assuming 
that  the  equilibrium  geometry  for  the  dimer  corresponds  to  a linear  con- 
figuration that  has  an  HF  hydrogen-bond  energy  of  6.  7 kcal/mole,  as 

34 

suggested  by  Coulson,  and  by  use  of  the  F-F  distance  in  the  HF  dimer 

the  value  2.  55  A predicted  from  an  ab  initio  self-consistent  field  (SCF)  study 

35 

carried  out  by  Del  Bene  and  Pople.  Their  value  of  5 reproduces  the 

known  experimental  dipole  moment  of  HF  at  the  equilibrium  separation. 

36 

Yarkony  et  al.  calculated  the  interaction  potential  between  two 

rigid  HF  molecules  by  means  of  the  Hartree-Fock  approximation,  Ab  initio 

35-39 

SCF  computational  studies  have  indicated  that  the  single-configuration 

SCF  theory  provides  a qualitative  correct  description  of  hydrogen  bonding. 

The  uncertainty  in  the  HF  dimer  bond  energy  makes  it  impossible  to  deter- 
mine if  the  SCF  theory  provides  a quantitative  description  of  hydrogen 
36  36 

bonding.  Yarkony  et  al.  published  total  energies  for  planar  and  non- 
planar  HF-HF  configurations.  The  coordinate  system  chosen  by  them  t<" 
describe  two  rigid  HF  molecules  with  an  internuclear  separation  of  1.723  bohr 
is  shown  in  Fig.  1.  This  coordinate  system  was  used  in  this  study  to  make 
it  much  easier  to  compare  parts  of  this  surface  with  those  of  Yarkony  et  al. 
for  R = 5,  6,  and  8 bohr.  The  resulting  contour  maps  are  shown  in  Figs.  2 
through  4.  For  the  planar  surface  calculations  at  R = 5,  6,  and  8 bohr,  the 


36 


least-square  deviations  a,  from  the  results  reported  by  Yarkony 
36 

et  al.,  are  as  follows;  R = 5aQ,  a = ± 2.  1 kcal/mole;  R = 6aQ,  o = ± 1.  3 
kcal/mole;  and  R = 8 aQ,  a = ± 0.25  kcal/mole.  This  semiempirical  surface 


compares  favorably  with  the  ab  initio  SCF  calculation  by  Yarkony  et  al. 


36 


Fig.  1.  Coordinate  system  used  for  HF-HF.  R is  the  distance  between 
the  two  HF  centers  of  mass.  The  angle  0 is  defined  with  re- 
spect to  the  indicated  z axis . <p  is  the  angle  of  rotation  of  the 
right-hand  HF  molecule  about  the  y-axis. 
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Fig.  2.  Contour  map  for  part  of  HF-HF  potential  surface  with  R = 5 ac 
and  0 = 0*  (planar).  8j  and  6 2 are  defined  in  Fig.  1.  The  con 
tours  are  in  units  of  kcal/mole,  relative  to  two  infinitely  sepa 
rated  HF  molecules.  Each  HF  molecule  has  a fixed  internu- 
clear  separation  of  1.737  a . 


Fig.  4.  HF-HF  potential  surface,  R = 8 aQ,  0 = 0*  (planar). 
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This  HF-HF  surface  and  the  SCF  surface  were  found  to  be  quite  repulsive 

at  values  of  R less  than  4.  5 bohr.  The  classical  trajectory  results  of  this 

study  indicate  that  part  of  the  surface  was  not  sampled  and  thus  is  not 

important  in  describing  either  vibrational  or  rotational  energy  transfer.  It 

is  indeed  very  encouraging  to  find  that  the  empirical  LEPS  parameters  used 

to  describe  the  energy  surfaces  for  HFg  and  FI^  systems  are  transferable 

to  the  four-atom  HF-HF  system.  Note,  however,  that  the  dipole-dipole 

term  played  a very  important  role  in  constructing  this  HF-HF  surface. 

36 

Yarkony  et  al.  showed  that  this  term  is  a necessary  part  of  any  reasonable 
analytic  fit  to  the  HF-HF  potential  energy  surface. 

The  values  of  R,  6 j,  9^,  and  <p  were  calculated  to  give  minimum  total 
energy  to  determine  the  dimer  geometry  for  the  surface.  This  surface  pre- 
dicts a planar  equilibrium  structure  with  0 j = 0*  , = 0°,  R(F-F)  = 2.9  A, 

36 

and  the  dimerization  energy  equal  to  2.7  ± 2 kcal/mole.  Yarkony  et  al. 

reported  0 j = 5%  0^  = 52®,  <f>  = 0®,  R(F-F)  = 2.  8 A,  and  the  dimerization 

energy  equal  to  4.5  kcal/mole.  The  reported  experimental  dimerization 

40-41 

energy  of  the  HF  dimer  is  (6  ± 1.5)  kcal/mole,  “ and  the  experimental 
equilibrium  geometry  is  0 j = 0°,  0^  = 60-70® , <p  = 0®,  and  R(F-F)  = 2.  8 A. 
The  results  of  all  these  studies  indicate  that  the  equilibrium  geometry  of  the 
HF  dimer  is  planar. 

43 

The  electron-gas  method  was  used  by  Parker  et  al.  to  calculate  the 
intermoiecular  potential  for  HF-HF  interactions.  In  this  method,  the  inter- 
molecular  energy  between  closed  shell  systems  is  the  sum  of  kinetic. 
Coulomb,  exchange,  and  correlation  energy  contributions.  The  detailed 
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method  for  calculating  intermolecular  potentials  between  closed-shell 

44  45 

systems  was  developed  by  Gaydaenko  and  Nikulin  and  Gordon  and  Kim 

and  was  modified  by  Rae^  and  Cohen  and  Pack.^  Parker  et  al.  ^ applied 

their  method  for  calculating  the  intermolecular  potentials  for  molecule- 

molecule  interactions  to  HF-HF  interactions.  The  results  of  this  study  are 

43 

compared  with  those  of  Parker  et  al.  and  with  the  accurate  SCF  calcula- 

36  43 

tions  of  Yarkony  et  al.  in  Figs.  5 through  7.  Parker  et  al.  added  the 

48 

leading  term  of  the  long-range  induction  (IND)  energy  to  the  SCF  part  of 
die  electron-gas  results  since  the  true  SCF  results  contain  induction  effects. 
The  method  described  here  to  calculate  the  intermolecular  potential  com- 
pares favorably  at  all  distances  and  angles  with  the  results  of  the  electron- 
gas  method  and  the  SCF  calculations.  The  present  method  and  the  electron- 
gas  plus  IND  method  do  not  contain  all  the  induction  and  charge  transfer 
effects  present  in  the  SCF  calculation  and,  therefore,  fail  to  describe 
accurately  die  hydrogen  bonding  region  shown  in  Fig.  5.  However,  this 
method  does  produce  enough  of  the  attractive  well  to  provide  a valid  calcula- 
tion for  energy  transfer  in  HF-HF  collisions. 
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HI.  CLASSICAL  TRAJECTORY  CALCULATIONS 


The  quasi-classical  procedure  described  by  Raff,  Thompson,  Sims, 

49 

and  Porter  was  used  to  examine  the  collision  dynamics  of  diatomic- 
diatomic  energy- transfer  processes.  The  Hamiltonian  that  describes  the 
four-particle  system  is  written  in  generalized  coordinates,  and  the  resulting 
set  of  18  Hamilton  equations  is  integrated  on  a CDC  7600  computer.  The 
initial  vibrational  states  of  HF(Vj,  J^)  and  HFtvg,  3^}  were  assigned  the 
vibrational  quantum  numbers  Vj  and  respectively.  Calculations  were 
made  with  values  of  v^  and  v ^ varying  from  0 through  6.  The  rotational 
quantum  numbers  Jj  and  J ^ were  Monte  Carlo  selected  for  the  calculations 
involving  v -*  v and  v-»  R energy-transfer  processes.  For  R -*  R,  T 
energy-transfer  processes,  the  rotational  quantum  number  Jj  was  assigned, 
and  the  rotational  quantum  number  J ^ was  Monte  Carlo  selected.  The  initial 
values  of  the  coordinates  in  each  collision  are  determined  by  means  of  a 
fixed  relative  translational  energy  of  the  reactants,  a chosen  initial  relative 
separation  between  the  centers  of  mass  of  the  diatomic  molecules;  chosen 
values  of  the  vibrational  quantum  numbers  Vj  and  v^;  and  randomly  selected 
values  of  the  rotational  quantum  numbers  J j and  J^,  the  impact  parameter, 
the  vibrational  phase  angles  of  both  molecules,  two  sets  of  orientational 
angles,  and  the  rotational  planes  of  both  moelcules.  Averaging  over  initial 
internal  states  of  HF(vj,  Jj)  and  HF(v2»  Jg)  was  carried  out  by  the  technique 
described  by  Porter,  Raff,  and  Miller.  ^ For  excited  vibrational  states, 
the  calculation  of  the  initial  values  of  the  momenta  P^(i  = 1, ...  ,9)  is  more 
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difficult  since  the  initial  amplitudes  of  both  molecules  are  no  longer  the 

classical  turning  points  and  the  internal  momentum  vectors  are  no  longer 

perpendicular  to  the  bond  axis.  For  this  case,  a method  based  on  Euler 

angles  was  used  and  that  was  adequately  described  by  Raff,  Thompson,  Sims, 
49 

and  Porter.  The  distance  R between  the  center  of  mass  of  the  two 
HF  molecules  was  taken  to  be  8 a.  u.  This  value  is  large  enough  to  ensure 
a negligible  initial  interaction  energy.  Relative  translational  energies  were 
assigned  values  that  ranged  from  0.  5 to  6 kcal/mole.  Approximately  200 
trajectories  were  run  for  each  initial  set  of  parameters.  A uniform  dis- 
tribution was  used  for  the  square  of  tile  impact  parameter  b.  All  runs  were 

made  with  a maximum  impact  parameter  of  2.  5 A.  The  step  size  was 
- 17 

5.7  x 10  sec.  The  integration  technique,  which  has  not  been  reported  in 

51  52 

the  literature,  proved  to  be  faster  than  the  Runge-Kutta-Gill  procedure, 

with  the  accuracy  of  the  integration  tested  by  changes  in  the  step  size  and  by 
integration  backwards  along  selected  trajectories.  As  an  additional  verifica- 
tion of  numerical  accuracy,  each  trajectory  was  checked  at  each  point  along 
the  trajectory  for  conservation  of  total  energy  and  total  angular  momentum. 
Computation  time  for  a single  trajectory  was  dependent  on  the  initial 
parameters  but  was,  on  the  average,  about  9 sec  per  trajectory  on  the  CDC 
7600  computer. 

The  final  properties  of  each  trajectory  were  analyzed  to  determine  the 
nature  of  the  collision,  i.  e. , tile  total  angular  momentum  and  vibrational- 
rotational  energy  of  both  HF  molecules.  The  partition  of  the  internal  energy 
between  vibrational  and  rotational  energy  was  determined  from  the  internal 
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energy  of  each  HF  molecule  and  its  total  angular  momentum.  The  actual 

technique  for  calculating  the  partitioning  of  the  internal  energy  has  been 

53 

described  previously  and  is  not  repeated  here.  The  v-»  v energy- 

transfer  cross  sections  were  calculated  by  Method  2 described  in  a previous 
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paper  by  Wilkins.  In  Method  2,  it  is  assumed  that  vibrational  energy  is 

transferred  in  every  collision  and  that  only  one  quantum  state  is  accessible 

in  flie  transfer  process.  The  energy- transfer  cross  sections  and  specific 

rate  coefficients  were  calculated  by  means  of  the  equations  given  by  Karplus, 
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Porter,  and  Sharma. 


IV.  RESULTS  AND  DISCUSSIONS 


A.  VIBRATIONAL -TO- ROTATIONAL  ENERGY  TRANSFER 


The  question  of  which  molecular  forces  cause  vibrational  energy 
exchange  in  HF-HF  collisions,  in  general,  has  not  yet  been  answered.  In 
hydrogen  halide  mixtures,  Chen  and  Moore17  demonstrated  in  HC1-HC1  and 
DC1-DC1  collisions  the  transfer  of  vibrational  energy  into  rotational  energy 
of  the  vibrationally  excited  hydrogen  halide  molecules.  In  analyzing  their 
relaxation  data  on  HC1  and  DC1,  Chen  and  Moore17  speculated  on  the  impor- 
tance of  such  strong  interactions  as  the  deep  attractive  well  in  the  potential 
energy  surface  caused  by  hydrogen  bonding.  The  formation  of  HF  dimers 
has  been  proposed  by  several  authors3-®*  1®* 11  to  explain  the  vibrational 
relaxation  of  HF(v  = 1)  by  HF(v  = 0).  The  dynamics  for  27,  000  HF-HF  col- 
lisions were  examined  for  the  potential  energy  surface  described  in  Section  II. 
A typical  HF-HF  bond  plot  is  displayed  in  Fig.  8.  A bond  plot  for  a four- 
atom  system  is  a plot  of  the  six  interparticle  distances  as  a function  of  time. 
The  motion  of  the  centers  of  mass  of  the  four -particle  system  as  it  moves  on 
the  potential  energy  surface  is  shown  in  Fig.  9.  The  typical  HF-HF  collision 
(Fig.  7)  indicates  that  the  two  F atoms  approach  one  another  until  repulsion 
reverses  their  course.  If  a collision  complex  were  formed,  then  R^jjcD^ 
in  Fig.  9 would  oscillate  back  and  forth,  which  provides  additional  evidence 
that  a complex  is  not'formed.  Since  there  are  six  degrees  of  freedom  in- 
volved in  constructing  the  potential  energy  surface,  it  is  not  possible  to  show 
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Fig.  8.  Bond  plot  of  HF(v  = 2)  + HF(v  = 0)  collisions.  R.. 

is  the  distance  between  the  ith  and  the  jth  par-  ^ 

tide,  t = 5.  7 X 10“  *7  sec. 
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Fig.  9*  Bond  plot  of  centers  of  mass  of  four -particle  system.  The  sym- 
bols A,  B,  C,  and  D are  defined  in  Fig.  1.  tQ  = 5.7  X 10*^  sec. 

I 


} 


the  nature  of  the  collision  on  a single  contour  plot  representing  only  two 
degrees  of  freedom.  This  trajectory  calculation  indicates  that  a collision 
complex  is  not  formed  for  the  typical  HF-HF  collision  in  the  temperature 
range  at  or  above  300  K.  The  present  calculations  indicate  that  the  HF-HF 
configuration  in  which  a hydrogen  bond  is  formed  is  not  a very  favorable 
configuration,  particularly  when  both  reagent  HF  molecules  are  rotating. 

The  trajectory  calculations  show  that  in  HF*-HF  collisions,  the  vibrational 

* 

energy  of  the  vibrationally  excited  HF  molecule  (HF  ) is  transferred  into 
rotational  energy  of  the  same  HF  molecule.  Specific  v -»  R rate  coefficients 
for  HF-HF  collisions  at  300  K are  given  in  Tables  II  through  IV.  The  qualita- 
tive results  of  these  v -•  R processes  are  as  follows: 

1.  The  vibrational  energy  of  the  vibrationally  excited  incident  HF 
molecule  is  transferred  into  rotational  energy  of  the  same  HF 
molecule.  There  is  practically  no  change  in  the  internal  energy 
state  of  the  target  HF  molecule. 

2.  In  these  v -*  R processes,  very  high  rotational  states  of  the 
incident  HF  molecules  are  populated  with  much  smaller  energy 
defects  than  would  have  been  predicted  if  both  reagent  and  prod- 
uct HF  species  were  assumed  rotationless. 

3.  The  v -*  R processes  occur  by  multiquantum  v transitions. 

HF(v  = 0)  species  can  be  formed  in  very  high  rotational  states. 

A loss  of  three  quanta  of  vibrational  energy  from  HF*  is  sufficient  to  form  an 
HF(v  = 0,  J = 24)  molecule,  with  an  energy  defect  of  125  cm  * (Table  IV). 


Table  II.  Vibration-to- Rotation  Detailed  Rate  Coefficients  for 
HF  (v.  = 1,  J.  = 2)  + HF  (v,  = 0)  at  T = 300  K 


by  (2  + 1)  e 


Table  III.  Vibration-to-Rotation  Detailed  Rate  Coefficients 
for  HF(vj  = 2,  Jj  = 2)  + HF  (v?  = 0)  Collisions 


'ill] 


1328. 0 
728.  0 


161.  0 


-370.  0 
-667.  0 
-1328. 0 
-1752.0 


1197.0 


406.  0 


-352.0 


-1075.0 


-1770.0 


-2427.  0 


2.  8 ± 1.  1 


2.  8 ± 1.1 


1.  4 ± 0.8 


2.  3 ± 1.0 


1.4  ± 0.8 


4.  7±  1.5 
3.  3 ± 1.  2 

2.  3 ± 

3.  3 ± 1.2 


2.  3 ± 


X = (vj,  Jj,  v2,  Vj,  Jj,  v2>  J2).  k^  is  the  total  rate  coefficient  for 

12  3 - 1 - 1 

v-»  R relaxation  in  units  of  10  cm  mole  sec.  To  obtain  rate  coef 

ficient  for  v-*  R relaxation  of  HF(v<  *2.  Ji  * 2)  by  HF(v,  = 0,  J_), 

-EV2=0,J2/kT  * 2 

multiply  kx  by  (2J?  + 1)  e /QR(v2»  J2^  where  QR(V2» 

is  the  rotational  partition  function. 
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Table  IV.  Vibration -to -Rotation  Detailed  Rate  Coefficients 
for  HF  (v.  = 3,  J4  = 2)  + HF  (v2  = 0) 

Collisions  at  T = 300  K 


V1 

J1 

v2 

J2 

B 

B 

/ 

V2 

J2 

cm 

k a 

K\ 

3 

2 

D 

J2 

2 

15 

0 

J2 

722.0 

2.  3 ± 1.0 

3 

2 

J2 

2 

14 

0 

J2 

178.0 

1.4  ± 0.8 

3 

2 

B 

J2 

2 

13 

0 

J2 

-333.  0 

2.  7 ± 1.1 

3 

2 

B 

3Z 

2 

12 

0 

J2 

-811.0 

3.2  ± 1.2 

3 

2 

3Z 

2 

11 

0 

J2 

-1254.0 

2.  7 ± 1.1 

3 

2 

B 

J2 

2 

10 

0 

J2 

-1662. 0 

6.  3 ± 1.6 

3 

2 

0 

J2 

1 

20 

0 

3Z 

437.  0 

3.  2 ± 1. 2 

3 

2 

0 

J2 

1 

19 

0 

J2 

-291.0 

0.9  ± 0.6 

3 

2 

0 

J2 

1 

18 

0 

J2 

-989.0 

1.8  ±0.9 

3 

2 

H 

J2 

1 

17 

0 

J2 

-1655.0 

2.  7 ± 1.  1 

3 

2 

B 

J2 

1 

16 

0 

J2 

-2289.0 

2.  7 ± 1.1 

3 

2 

0 

J2 

0 

26 

0 

J2 

1962.  0 

1.4  ± 0.8 

3 

2 

0 

J2 

0 

25 

0 

h 

1030.0 

4.  5 ± 1.4 

3 

2 

0 

J2 

0 

24 

0 

h 

125.  0 

2.  7 ± 1.  1 

3 

2 

0 

3Z 

0 

23 

0 

3 2 

-752. 0 

2.  3 ± 1.0 

3 

2 

0 

J2 

0 

22 

0 

J2 

-1602.0 

2.3  ± 1.0 

3 

2 

0 

J2 

0 

21 

0 

J2 

-2421.0 

1.8  ± 0.9 

= (Vj,  v2»  v(»  v2*  tot**  rate  coefficient  for  v -*  R relaxation  is 

in  units  of  10*^  cm^  mole"*  sec"*.  To  obtain  rate  coefficient  for 
v -*  R relaxation  of  HF(Vj  = 3,  Jj  = 2)  by  HF(v2  = 0,  J 2),  multiply  by 

"Ev  = 0,  J /kT 

(2J2  + 1)  e /Qr(v2»  J2)»  where  Qj^(v2»  Jj)  **  t*le  rotationa^ 

partition  function. 
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The  ( v ' , j')  states  with  the  smallest  energy  defects  to  be  populated  by  v -•  R 
processes  in  HF-HF  collisions  are  listed  in  Tables  V and  VI.  The  v -*  R pro- 
cesses are  quite  complex  (Tables  V and  VI).  When  an  HF(v  = 3)  molecule 
collides  with  an  HF(v  = 0)  molecule,  the  HF(v  = 3)  molecule  can  transfer  one, 
two,  or  three  quanta  of  vibrational  energy  into  rotational  energy  of  the  same 
HF  molecule  by  an  intramolecular  energy-transfer  mechanism  (Table  V). 
When  both  HF  molecules  are  vibrationally  excited,  either  or  both  HF*  mole- 
cules are  capable  of  converting  one  or  more  quanta  of  vibrational  energy  of 
either  or  both  HF*  molecules  by  intramolecular  energy  transfer  (Table  VI). 

It  is  probable  that  v - R processes  are  required  to  explain  the  high  J states 
observed  in  the  pure  rotation  laser  described  by  Deutsch.2® 

This  is  the  first  time,  to  our  knowledge,  that  v - R processes  have 
been  predicted  by  a three-dimensional  trajectory  calculation  to  be  very  im- 
portant in  HF-HF  collisions.  It  is  not  the  first  time,  however,  that  the  role 
of  v **  R mechanisms  in  atom -diatomic  collisions  or  diatomic -diatomic  colli- 
sions have  been  recognized.  Kelley^  investigated  the  role  of  rotational 
motion  in  the  de-excitation  of  vibrationally  excited  diatomic  molecules  with 
an  atom.  Application  of  his  impulsive  collision  theory  indicates  that  intra- 
molecular v -*  R processes  should  be  very  important  in  HF-HF  collisions. 

57 

Dillon  and  Stephenson  reported  that  multiquantum  rotational  transitions 

play  an  important  role  in  vibrational  energy  exchange  involving  diatomic - 

diatomic  collisions  by  allowing  large  vibrational  energy  defects  to  be  ab- 

58 

sorbed  by  the  rotational  degrees  of  freedom.  McGuire  and  Toennies  used 
the  couple-state  approximation  to  calculate  vibrational  de-excitation  cross 
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Table  VI.  Most  Probable  (v/,J/)  States  to  be  Populated  by 
Vibrational-to-Rotational  Energy  Transfer  in 
HF(vj)  + HF  (v^)  Collisions  (Continued) 


vi 

v2 

B 

B 

v. 

v2 

B 

5 

5 

4 

14 

4 

14 

5 

5 

3 

20 

4 

14 

5 

5 

3 

20 

3 

20 

5 

5 

2 

24 

4 

14 

5 

5 

2 

24 

3 

20 

5 

5 

2 

24 

2 

24 

5 

5 

1 

28 

4 

14 

5 

5 

1 

28 

3 

20 

5 

5 

1 

28 

2 

24 

5 

5 

1 

28 

1 

28 

5 

5 

0 

32 

4 

14 

5 

5 

0 

32 

3 

20 

5 

5 

0 

32 

2 

24 

5 

5 

0 

32 

i 

28 

5 

5 

0 

32 

0 

32 

sections  for  the  He-H ^ system.  They  found  that  v -*  R energy  transfer  plays 
a major  role  in  relaxation  by  Ar.  The  role  of  v -*  R transitions  has  been 
discussed  in  connection  with  experimental  observations  for  the  He-^ 
system^”***  and  He-D systems.  Nazar  et  al.  ^ postulated  that  highly 

rotationally  excited  hydrides  can  relax  by  cascading.  This  process  involves 
rapid  vibrational  exchange  in  which  the  vibrationally  excited  hydride  (the 
donor  molecule)  is  de -excited  by  a single  quantum  of  vibrational  energy  with 
retention  of  its  initial  rotational  quantum  number.  Although  the  Nazar  et  al. 
observation  does  not  concur  with  our  v -*  R prediction  that  a large  change  in 
AJ  occurs  by  an  intramolecular  energy  mechanism  in  the  donor  molecule, 
their  experiment  indicates  that  this  rapid  energy  exchange  is  indicative  of  a 
resonant  v -*  R energy-transfer  process.  The  theory  of  Dillon  and  Stephenson, 
however,  does  concur  with  finding  of  this  study  that  there  are  efficient  v -•  R 
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processes  in  which  large  vibrational  energy  defects  are  absorbed  by  multi- 
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quantum  changes  in  the  rotational  degree  of  freedom.  Streit  and  Johnston 
in  their  work  on  the  quenching  of  vibrationally  excited  OH  by  several  chap- 
erones suggested  that  in  such  collision  processes  a considerable  amount  of 

vibrational  energy  is  absorbed  into  rotation  in  order  to  reduce  the  energy 

65 

defect.  This  idea,  according  to  Streit  and  Johnston,  requires  a theory  in 

which  A J transitions  are  responsible  for  the  very  efficient  vibrational  relaxa- 

67 

tion  of  HF(v  < 9).  Downey  et  al.  found  intramolecular  v -*  R energy-transfer 
mechanisms  to  be  important  in  fee  relaxation  of  vibrationally  excited  OH  mole- 
cules by  Ar  atoms.  This  intramolecular  v -*  R transfer,  according  to 
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Downey  et  al.  f**  must  be  the  dominant  channel  for  relaxation  in  experiments 
involving  a light  rotor  and  Ar,  or  a light  rotor  and  Krypton. 

In  Tables  VII  through  IX  are  listed  detailed  rate  coefficients  at  T = 300  K 
for  the  v -*  R energy -transfer  processes  HF(Vj)  + HFfv^)  ■*  HF(v')  + HF(v^). 

, / / are  easily  converted  to  probabilities 

vrWvz  13  , 

/ /by  dividing  k / / by  the  gaskinetic  rate  9 X 10  cm 

vrv2;vrv2 


vrv2:vrv2 


The  v -*  R rate  coefficients  k 
P 

mole ' * sec*!  at  room  temperature.  There  is  a probability  of  0.2  for  most 
single -quantum  transitions  and  a probability  of  0.  i for  most  multiquantum 
transitions  from  any  vibrational  level  v (Table  VII).  The  overall  rate  coeffi- 
cients tabulated  in  Tables  VII  through  IX  indicate  a slight  increase  with  in- 
creasing vibrational  quantum  number  v.  These  overall  rate  coefficients  are 
not  comparable  to  the  effective  rate  coefficients  measured  in  the  laboratory. 
The  rate  coefficients  listed  in  Tables  VII  through  IX  are  for  all  v -*  R pro- 
cesses from  a given  v^  level  to  a final  Vj  level,  regardless  of  the  final  Jj 
state.  The  rate  coefficients  for  v -*  R transfer  processes  in  which  the  final 
J*  states  are  given  in  Tables  II  through  IV.  When  the  detailed  v-»  R rate 
coefficients  listed  in  Tables  II  through  IV  are  used  in  a rotational 
nonequilibrium  model  along  with  R — R* , T rate  coefficients  presented  here, 
temperature-dependent  rate  coefficients  are  calculated  for  HF(v  = 1)  + 

HF(v  = 0)  collisions  that  concur  with  measured  temperature -dependent  effec- 
tive v — R rate  coefficients.  The  actual  calculations  will  be  described  by 
24 

Kwok  and  Wilkins  at  a later  time.  There  is  an  increase  in  the  number  of 
v *•  R channels,  and  a decrease  in  the  specific  rate  coefficients  for  v -•  R 
processes  with  increasing  v (Tables  VIII  and  IX).  The  multiquantum  v -*  R 


- 51  - 


I 


Table  VIII. 


i 


Rate  Coefficients  *or  Vibrational -to -Rotational 

Energy  Transfer  in  HF^Vj)  + HF(v2)  Collisions  at  T = 300  K 


V1 

v2 

V2 

k / / 

vi-v2-vrv2 

k a 
vR 

1 

1 

i 

4.  1 ± 1.3 

1 

1 

0 

1.5  ± 0.9 

5.  6 ± 0.3 

2 

2 

2 

i 

2.  8 ± 1.1 

2 

2 

2 

0 

2.0±  1.0 

2 

2 

1 

1 

0.7  ± 0.2 

2 

2 

1 

0 

0.  5 ± 0.  1 

2 

2 

0 

0 

1.  0 ± 0.2 

7.  0 ± 0.3 

3 

3 

3 

2 

1.7  ± 0.2 

3 

3 

3 

1 

1.  5 ± 0.2 

3 

3 

3 

1.5  ± 0.2 

3 

3 

2 

0.7  ± 0.  2 

3 

3 

2 

0.  5 ± 0.  1 

3 

3 

2 

0.3  ± 0.  1 

3 

3 

1 

1 

0.  1 ± 0.  1 

3 

3 

1 

1 ; 1 

0.6  ± 0.  1 

3 

3 

0 

0.  6 ± 0.  1 

7. 4 ± 0.3 

4 

4 

3 

4 

1.7  ± 0.2 

4 

4 

2 

4 

0.  9 ± 0.2 

4 

4 

1 

4 

0.7  ± 0.2 

4 

4 

4 

1 . 1 ± 0.2 

4 

4 

3 

3 

0.2  ± 0.  1 

4 

4 

2 

3 

0.6  ± 0.  1 

4 

4 

1 

3 

0.4  ± 0.  1 

4 

4 

0 

3 

0.2  ± 0.  1 

4 

4 

2 

2 

0.  1 ± 0.  1 

4 

4 

1 

2 

0.  2 ± 0.  1 

4 

4 

0 

2 

0.  5 ± 0.1 
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Table  VIII.  Rate  Coefficients  ky  v2,  vl>  v2  for  Vibrational-to-Rotational 
Energy  Transfer  in  HF(vj)  + HF(v2)  Collisions  at  T = 300K 


(Continued) 


c / / is  in  units  of  10*3  cm3  mole"  sec' 

vl,v2;vl’v2 
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Table  IX.  Rate  Coefficient  k_r , / • for  Vibrational-to -Rotational 

M,  v2;  V4,  v2 

Energy  Transfer  in  HF(vj)  + HF(v2)  Collisions  at  T = 300  K 


v2 

v2 

k / / 

vl*  V2*  vi’  V2 

k * 
vR 

1 

i . 

1 

1.6  ± 0.2 

1 

0 

1 

1.7  ± 0.2 

1 

2 

0 

2.  0 ± 0.2 

1 

0 

0 

1.  1 ± 0.2 

6.  4 ± 0.  3 

1 

2 

1 

1.2  ± 0.2 

1 

1 

1 

0.  5 ± 0.  2 

1 

0 

1 

1.4  ± 0.2 

1 

3 

0 

1.8  ± 0.2 

1 

2 

0 

0.7  ± 0.2 

1 

0 

0 

0.8  ± 0.2 

6. 5 ± 0. 3 

1 

3 

1 

1.0±  0.2 

1 

2 

1 

0.6  ± 0.2 

1 

1 

1 

0.4  ± 0.  1 

1 

0 

i 

1.3  ± 0.2 

1 

4 

0 

1.5  ± 0.2 

1 

3 

0 

0.8  ± 0.2 

1 

2 

1 1 

0.6  ± 0.2 

1 

0 

0 

■ W 'era 

6.7  ±0.3 

1 

4 

1 

; ' ' . ' ; / 

1 

3 

1 

1 

2 

1 

0.3  ± 0.  1 

1 

1 

1 

0.  3 ± 0.  1 

1 

0 

1 

1.0  ± 0.2 

1 

5 

0 

1.4  ± 0.2 

1 

4 

0 

0.8  ± 0.2 

1 

3 

0 

0.  5 ± 0.  1 

i 

2 

0 

1 

0 

0 

■BIB 

6.  9 ± 0.3 

*k  / / is  in  units  of  10^  cm^  mole"*  sec"*. 

VT  v2;v  1'  v2 
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energy-transfer  processes  are  slightly  less  probable  than  the  single -quantum 

17 

processes.  Chen  and  Moore's  model  of  their  measurement  for  vibrational 
relaxation  of  HC1  by  HC1  and  DC1  by  DC1  support  our  prediction  that  rota- 
tional up-pumping  occurs  as  the  vibrational  energy  is  reduced  by  deactiva- 
tion in  HF-HF  collisions. 


B. 


VIBRATIONAL -TO- VIBRATIONAL  ENERGY  TRANSFER 


The  v -*  v energy-transfer  processes  are  one  of  the  kinds  of  energy- 
transfer  processes  that  are  responsible  for  the  excitation  and  de -excitation 
of  vibrationally  excited  HF  species  in  HF  chemical  lasers.  It  is  important 
in  evaluating  the  performance  of  HF  chemical  lasers  to  have  temperature- 
dependent  rate  coefficients  for  the  v -*  v energy -transfer  processes.  Since 
the  energy  mismatch  AE  is  small  for  these  collisions,  there  has  been  specu- 
lation that  the  energy-transfer  rates  for  the  v -*  v processes  are  very  fast. 

Rate  coefficients  for  room  temperature  have  been  provided  by  several  experi- 

1 6 - 8 

mental  studies  * for  endothermic  processes  (reading  from  left  to  right) 


^1 

HF(Vj,  Jj)  + HF(v2  =0,  J2)  g=t  HF(vJ  = Vj 


i,j;)  + HF(v'  = 1,J')  (1) 


with  Vj  = 2 through  5.  If  it  is  assumed,  as  it  commonly  is,  that  the  net 
change  in  rotational  energy  across  the  endothermic  processes  is  zero  and, 
if  equilibrium  constants  are  used  to  calculate  the  rate  coefficients  in  the  exo- 
thermic direction,  then  rate  coefficients  for  some  of  these  processes 
as  high  as  10  Torr  sec"  at  room  temperature  are  obtained. 
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Another  problem  associated  with  the  experimental  measurements  is 

* 

that  in  the  measurements  it  has  not  been  possible  to  separate  the  contribu- 

- tions  from  the  v -*  v and  v -*  R,T  energy-transfer  processes.  In  addition  to 

this  problem,  there  is  the  problem  of  not  knowing  experimentally  the  actual 
energy  defects  AE  for  the  endothermic  processes  given  in  Eq.  (1).  The 
main  advantage  of  this  trajectory  study  is  that  the  rate  coefficients  for  both 
the  endothermic  and  exothermic  v -»  v energy-transfer  processes  can  be 
calculated.  Rate  coefficients  for  room  temperature  v -*  v energy^transfer 
processes  obtained  from  this  trajectory  study  are  given  in  Table  X.  In 
column  6,  of  Table  X,  is  given  the  energy  mismatch  AE  obtained  from  the 
relationship  between  the  rate  coefficients  for  the  endothermic  and  exothermic 
processes,  i.e..,  k } = kj  exp(AE/RT).  The  trend  of  AE  is  not  apparent 
from  this  calculation  since  there  is  a 25%  uncertainty  in  these  calculated 
rate  coefficients.  However,  it  is  clear  that  the  v -•  v processes  are  very 
much  near  resonant.  In  column  7 are  listed  the  energy  defects  obtained 
under  the  assumption  that  both  reagent  and  product  HF  species  have  zero 
rotational  energy.  The  large  rate  coefficients  predicted  experimentally  for 
the  exothermic  processes  can  be  explained  only  by  application  of  these  large 
energy  defects.  The  trajectory  calculations  indicate  that  such  large  energy 
defects  are  not  correct  and  that  the  rotational  states  of  reagent  and  product 
HF  species  must  be  considered  to  calculate  a realistic  energy  defect  AE. 

On  the  basis  of  the  energy  defects  tabulated  in  column  7,  the  rate  coefficients 
for  v -*  v energy-transfer  processes  would  be  expected  to  decrease  with  in- 
creasing  v.  The  trajectory  calculations  predict  (column  5)  that  the  rate 
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Table  X.  Rate  Coefficients  for  Vibrational-to -Vibrational  Energy  Transfer 
HF(v.)  + HF(v_)  -*  HF(v.  - 1)  + HF(v_  + i) 


<4 


w e 
<i  5 

m on  M*  <n  cm 

(s-  if  O sO  (M 

»*  Bi  in  vO  co 

in  m ^ id  (M 
r-  ^ o vo  <m 
^ co  m so  oo 

i i i i i 

vO  no  so  in  in 

4H  ^ 

1 1 1 1 1 

168 

AE,  b 

-1 

cm 

N f-  00  **  ■" 

n-  oo  co  oo  oo 

CM  f"  00  "■*  "< 
l>  CO  00  CO  CO 

1 1 1 1 1 

m 

i 

+13 

▼4 

u 

+ ® 
►"  7 

O'lfNO'N 

7 00  af*  CM  O' 

oo  cm  cm 

- ® 

- 1 

-H  *H  -H  -H  -H 

0 ±2 
± 1. 
± 1. 
± 1. 
± 0. 

4H  4H  ^ O O 

-H  +»  -H  -H  -H 

• 

-H 

V2SV 

3 

m n 

m cm  an  O'  cm 

oo  ici  m " " 

• OOOOOI' 

CM  • • 

r-  m cm  "< 

'OO'IBIMN 

vO  IB  N "0 

CM 

• 

vO 

. u 

> ■"« 

4*  o 

+ 

(M 

> 

1 

1 

1 

1 

1 

N IB  4*  m nO 

m m*  m so  r- 

CM 

> 

■h  n cn  4*  m 

o o o o o 

•«-<  cm  on  M*  m 

IM 

CM 

> 

o o o o o 

•HNiBfifl 

NlB^lB'O 

> 

cm  m ■»(*  »n  vO 

1 

1 

1 

1 

1 

im  cn  M*  m so 

m 

V 

<4 

h 

C 

•H 

l4 

-M 

43 

O 


c 

V 

60 

a 

v 


o 

m 


in 

o 

o 

m 


>S 

43 

■O 

V 


I 

<u 

43 


c ® 

* M 
4)  3 

H -M 

43  i4 
<6  h 

*1 

c +* 

•H  M 
V 4) 

•3$ 

<4  O 
> >, 

£9 

H +j 
(4 

• -•-> 
Sd  C 
^ 4) 
O •-* 
O O 

CO  £ 
*«-l 

II  4) 

H 8 

<4  r 


aa 

C 

0 


O bo 
M U 

V V 
N G 

« ® 

•°-a 

o c 
+»  o 
x)  *> 
41  (4 

es 

3 M 

S 2 
- 5 

m N 

•H 

g > 

O (4 

•H  !2 

So 

4) 

* *0 
_ 4) 

: & 
O 3 
h « 

u aa 
<4  (4 

tK  4) 

60  M 
M (4 

V _ 
C • 

® .a 


5 ■a 


(4 

3 

u 

3 

u 

>s 

0 

u 

4) 

•n 

<4 

H 


u 

4) 

a 

aa 


- 58" 


hi 


coefficients  decrease  with  v for  both  the  endothermic  and  the  exothermic 
v -*  v energy-transfer  processes.  From  this  trajectory  calculation,  it  is 
predicted  that  only  single -quantum  processes  are  important  in  the  v -»  v 
intermolecular  energy -transfer  processes  and  that  the  v dependence  in  the 
endothermic  direction  is  approximately  given  by  k n.  , . = 2 “ n . . 

Ld  p U 1 1 p 1 

for  v > 2.  In  the  exothermic  direction,  the  v dependence  is  k . , , = 

c v,v;v-l,v+l 

1 - v 

2 kj  j,Q  2 for  v > 1.  Trajectory  calculations  were  performed  for  T = 300> 

500,  1000,  and  1500  K.  Analysis  of  the  results  indicates  a temperature  de- 
0 5 

pendence  of  T * for  both  the  v -*  v and  v -»  R energy-transfer  processes. 

In  Table  XI  are  given  the  measured  (v  -*  v + v — R)  rate  coefficients, 
the  v -*  R rate  coefficients  for  an  assumed  v dependence,  the  reduced  endo- 
thermic v -•  v rate  coefficients,  v “♦  v rate  coefficients  obtained  from  this 
trajectory  study,  and  v -*  v rate  coefficients  calculated  by  Shin***  with  his 
one -dimensional  analytical  model.  In  columns  2 through  4 are  listed  the 
measured  rate  coefficients  for  the  (v  -*  v + v -*  R)  energy -transfer  processes. 

In  column  5 are  listed  the  v -*  R rate  coefficients  with  kv_>^  = v^'  ^ kY-*? 

v l-*0 

assumed.  The  2.  3 power  of  v was  recommended  by  Kwok  as  a result  of  his 
flowtube  studies  on  the  HF  + HF  system.  In  columns  6 through  8 are  listed 
the  v “•  v rate  coefficients  obtained  from  the  measured  (v  -*  v + v -»  R)  rate 
coefficients  by  subtracting  from  the  sum  the  value  given  in  column  5 for  the 
assumed  v -»  R rate  coefficient.  Although  the  measurements  for  the  (v  -*  v + 
v -»  R)  rate  coefficients  agree  within  a factor  of  two,  the  uncertainty  in  all  of 
the  measured  (v  -»  v + v -*  R)  rate  coefficients  is  too  large  (in  excess  of  25%) 
to  deduce  the  correct  v dependence  for  the  v -•  v energy-transfer  processes. 


1 1 
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If  the  v -»  v rate  decreases  with  increasing  v,  as  predicted  from  this  trajectory 
study  (column  9),  then  the  (v  -*  v + v -•  R)  measured  rates  at  the  high  v levels 
(v  greater  than  4)  must  reflect  only  v -•  R energy-transfer  processes.  The 
slow  v “•  v pumping  rates  and  the  fast  v -•  R de -excitation  rates  with  increas- 
ing v could  account  for  the  large  decrease  in  measured  excited  HF  species 
concentrations  in  the  higher  vibrational  levels  (v  greater  than  6).  The  v -•  v 
rate  coefficients  obtained  from  this  trajectory  study  (column  9)  can  be  com- 
pared with  those  reported  by  Shin*^  (column  10)  with  his  one -dimensional 
analytical  model.  The  v — v rates  reported  by  Shin^  indicate  a slight  in- 
crease from  v = 2 to  v = 4,  then  a decrease  from  v = 4 to  v = 5;  this  increase 
from  v = Z to  v = 4 is  not  seen  in  results  obtained  from  the  trajectory  study. 
The  arrested  rotational  model  used  by  Shin  in  which  the  mechanism  for  v -*  v 
energy  transfer  in  HF-HF  collisions  is  the  result  of  strong  hydrogen -bond 
attraction  was  not  observed  in  the  27,  000  trajectories  analyzed  in  this  study. 
These  trajectories  show  that  most  of  the  energy  mismatch,  corresponding 
to  rotationless  reagent  and  product  HF  molecules,  enters  rotational  motion 
of  the  colliding  HF  molecules.  This  study  indicates  that  dimer  formation  at 
temperatures  at  and  above  room  temperature  will  not  occur  if  both  reagent 
HF  molecules  have  initial  rotational  energy,  since  the  dipole -dipole  forces 
are  strongly  attractive  only  when  the  HF  molecules  are  oriented  in  an  unfavor- 
able configuration.  If  HF  dimers  were  to  form*  then  the  v *■'  v rates  would 
increase  with  increasing  v because  the  density  of  states  in  the  (HF^  transi- 
tion complex  would  increase  with  increasing  vibrational  excitation.  This 
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trajectory  study  and  the  reduced  v -»  v rate  coefficients  obtained  from  the 
measured  (v  -*  v + v -*  R)  rate  coefficients  do  not  seem  to  reflect  this 
increase. 

C.  ROTATIONAL- TO -ROTATIONAL,  TRANSLATIONAL 
ENERGY  TRANSFER 

68 

In  a recent  article,  Gur'ev  et  al.  obtained  rate  coefficients  for  rota- 
tional relaxation  of  HF(v  = 0,  J = 8)  by  several  chaperone  gases  by  measur- 
ing the  relaxation  losses  that  occur  when  light  is  made  to  pass  through  a 
resonantly  absorbing  medium.  Their  rate  coefficients  for  rotational 
relaxation  of  HF(v  = 0,  J = 8)  by  M,  where  M is  an  atom  or  molecule,  cor- 
respond to  collision -induced  relaxation  probabilities  of  0.  1 for  M = HF, 

3 X 10-3  for  M = H2  or  D^,  and  <5  X 10-^  for  M = He,  Ar,  Kr,  or  Xe. 

68 

Gur'ev  et  al.  point  out  that  such  low  probabilities  call  for  a more  thorough 

analysis  of  vibrational- rotational  relaxation  of  HF  molecules  in  comparison 

with  the  usual  assumption  that  R -*  R,  T relaxation  is  instantaneous. 

69 

Vasil'ev  et  al.  showed  that  resonance  transfer  of  rotational  energy  from 
HF  molecules  to  H ^ or  plays  an  important  role  in  rotational  deactivation 
of  rotationally  excited  HF  molecules  with  J 5 7.  They  found  that  the  rate  of 
rotational  deactivation  of  HF(v/,  Jr)  by  either  He  or  Xe  decreases  with  increas- 
ing j'  and  that  more  than  100  collisions  are  required  for  complete  deactiva- 
tion of  rotationally  excited  HF  molecules  formed  in  the  F + H^,  reaction. 
Vasil'ev  et  al.^  also  found  that  R -•  R,  i.e.  , resonance  transfer  in  HF  + 
and  HF  + systems  enhance  the  rotational  deactivation  rates  by  an  order 
of  magnitude,  compared  with  the  deactivation  by  He  or  Xe  atoms. 
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The  rate  coefficients  for  this  study  for  the  rotational  relaxation  of 
HF(Vj  = 0,  Jj  = 5,  10,  15)  by  HF(v£  = 0,  3^)  molecules  are  shown  in 
Figs.  10  through  12.  These  rate  coefficients  contain  contributions  from 
both  R -*  R and  R -*  T energy-transfer  processes.  The  R -*  T mechanism 
with  AJ  = -1  is  the  main  mechanism  for  de -excitation  of  rotationally  excited 
HF  species.  The  contribution  from  the  R -»  R mechanism  for  AJ  = -1  is 
approximately  one -third  of  the  total  contribution  to  the  (R  -*  R + R -*  T)  rate 
coefficients.  At  low  J,  de-excitation  of  rotationally  excited  HF  molecules 
by  multiquantum  rotational  transitions  is  much  more  probable  than  at  high 
J (Fig.  12), where  de-excitation  with  AJ  = -2  is  more  than  one  order  of  mag- 
nitude smaller  than  de-excitation  with  AJ  = -1.  The  total  rate  coefficients 
for  de -excitation  of  rotationally  excited  HF(Vj  = 0 and  1,  J^  = 5,  10,  15)  by 
HFOv^  = 0,  J molecules  are  shown  in  Figs.  13  and  14.  It  can  be  seen  from 
these  curves  that  the  rate  coefficients  for  de -excitation  from  the  higher  rota- 
tional levels  are  much  smaller  than  rate  coefficients  for  de -excitation  from 

the  lower  J levels.  These  results  concur  with  those  reported  by  Polanyi  and 
18  19 

Woodall  and  Ding  and  Polanyi,  who  showed  that  for  hydrogen  halides  the 
probability  of  rotational  deactivation  decreases  with  increasing  rotational 
quantum  number  J.  The  tendency  for  the  high  J levels  to  relax  less  rapidly 
than  the  lower  J levels  is  a consequent  of  the  large  spacings  between  adja- 
cent rotational  states  of  high  J.  Most  of  the  rotational  transitions  observed 
20 

by  Deutsch  had  J values  greater  than  15,  which  indicates  that  the  high 
rotational  states  had  not  Boltzmannized.  The  nonequilibrium  rotational 
distribution  occurs  because  the  rotational  level  spacings  become  so  large 
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Fig.  10.  Rate  coefficients  for  R -*  R,  T energy  transfer  of  HF(v.  = 0,  J.  = 5) 
by  HF(v2  = 0,  J2).  HF(v«  = 0,  Jj  = 5)  + HF(v,  = 0,  J2)  - 
HF(Vj  = 0,  J'2  =^x  + AJ)  + HF(vJ  = 0,  J'2). 
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Fig.  12.  Rate  coefficients  for  R **  R,  T energy  transfer  of  HF(vj  = 0,  Jj  = 15) 
by  HF(v2  = 0,  J2).  HF(V4  = 0,  Jj  = 15)  + HF(v,  = 0,  J2)  - 
HF<v'j  = 0,  Jj  + AJ)  + HF(v'2  = 0,  J'2). 
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Fig.  13.  Total  rate  coefficients  for  rotational  de-excitation  of 
HF(vj  = 0,  Jj  =5,  10,  and  15)  by  HF(v2  = 0,  J2). 
HF(vj  = 0;Jj)  +HF(v2  = 0,J2)-  HF(vj  = 0,Ji  + AJ) 

+ HF(v2  = 0,J'). 
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that  a few  collisions  no  longer  induce  rotational  relaxation  by  AJ  = -1 
transitions.  The  transition  probabilities  for  rotational  de-excitation  of 
HF(Vj  =0,  Jj)  by  HF(v2  = 0,  molecules  at  T = 300  K are  given  in 
Table  XII.  Rotational  relaxation  out  of  the  high  J levels  into  the  low  J 
levels  by  multiquantum  transitions  is  a very  inefficient  process.  A transi- 
tion from  J = 15  to  J = 12  occurs  about  one  in  a thousand  collisions.  A tran- 
sition from  J = 3 to  J = 2 occurs  about  one  in  three  collisions,  whereas  a 
transition  from  J = 3 to  J =0  occurs  about  one  in  fifty  collisions.  Multi - 
quantum  J transitions  are  much  more  probable  at  the  low  J levels  than  at 
the  high  J levels . Even  at  low  J levels , the  probability  of  multiquantum 

transitions  are  at  least  an  order  of  magnitude  smaller  than  the  probability 

70 

of  a single -quantum  A J transition.  Hinchen  and  Hobbs  measured  rota- 
tional relaxation  effects  in  HF  + M collisions  by  means  of  a double  resonance 
infrared  experiment.  Their  R -*  R,T  rates  at  low  J values  are  at  least  one 

order  of  magnitude  larger  than  the  values  predicted  in  this  study.  The  rate 

70 

model  used  by  Hinchen  and  Hobbs  to  calculate  the  specific  rate  coefficients 
did  not  include  v -*  R and  R -*  v processes  found  to  be  very  important  in 
HF-HF  collisions  in  this  study. 
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Table  XII.  Rotational  Transition  Probabilities  for 
HF-HF  Collisions  at  T = 300  K 


AJ  = -1 

AJ  = -2 

AJ  = -3 

0.45 

0.40 

0.  065 

0.  37 

0.  055 

0.  022 

0.35 

0.047 

0.020 

0.32 

0.040 

0.019 

0.30 

0.  030 

0.  017 

0.28 

0.025 

0.015 

0.25 

0.017 

0.014 

0.22 

0.  015 

0.012 

0.20 

0.010 

0.010 

0.  18 

0.009 

0.  008 

0.  15 

0.  008 

0.  007 

0.  14 

0.  007 

0.005 

0.  12 

0.006 

0.003 

0.  10 


0.  005 


0.001 


V.  SUMMARY 


The  v -*  v processes  occur  by  intermolecular  energy  exchange 
mechanisms  in  which  the  most  probable  paths  involve  single  vibrational 
quantum  transitions.  The  v -»  v rate  coefficients  decrease  with  increasing 
v,  and  most  v -»  v processes  are  near  resonant  with  an  energy  defect  of  less 
than  100  cm"^.  Both  the  vibrational  and  the  rotational  states  of  the  reagents 
and  product  species  must  be  used  in  calculating  energy  defects.  When  proper 
energy  defects  are  used,  both  endothermic  and  exothermic  HF-HF  v -*  v rate 
coefficients  decrease  with  increasing  v,  and  the  values  of  the  rate  coeffi- 
cients in  the  endothermic  or  the  exothermic  direction  will  not  exceed  the 
gas  kinetic  value. 

The  v -*  R processes  occur  by  intramolecular  energy-transfer  mecha- 
nisms. The  vibrationally  excited  HF  molecule  converts  one  or  more  quanta 
of  its  vibrational  energy  into  rotational  energy.  In  an  HF  system,  it  is  pre- 
dicted that  v -♦  v up-pumping  coupled  with  v -*  R de -excitation  with  large  AJ 
changes  provide  the  mechanisms  necessary  to  explain  the  high  J states  ob- 
served by  Deutsch  in  his  HF  rotational  laser.  The  v -*  R rate  coefficients 
increase  with  increasing  v.  The  vibrationally  excited  HF  rotor  takes  the 
initial  energy  mismatch  AE,  corresponding  to  rotationaless  HF  molecules, 
away  through  v -*  R nonresonant  energy-transfer  processes.  This  initial 
energy  defect  goes  into  rotational  motion  of  the  colliding  molecules.  For 
many  of  the  important  v -*  R energy -transfer  processes,  the  final  energy 
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defect  for  the  process  is  less  than  200  cm'l  The  multiquantum  v -»  R 
transitions  provide  more  ways  to  distribute  the  vibrational  energy  of  the 
vibrationally  excited  HF  molecules  into  rotational  energy,  i.e.,  into  very 
high  rotational  quantum  states. 

Many  of  the  R -»  v processes  are  fast  because  of  the  small  energy  defect 
(about  200  cm  *)  found  for  the  corresponding  v -*  R processes. 

Both  R -*  R and  R -*  T processes  occur  in  the  rotational  de-excitation  of 
rotationally  excited  HF  molecules.  The  probability  of  rotational  de-excitation 
decreases  with  increasing  rotational  quantum  numbers  J.  The  single  rota- 
tional quantum  transitions  are  more  probable  than  the  multiple  rotational 
quantum  transitions. 

HF  dimers  do  not  have  to  be  formed  at  room  temperature  to  explain  the 
fast  v -*  R self-relaxation  rates  measured  by  various  experimental  techniques. 
In  a subsequent  report,  the  temperature -dependent  rate  coefficients  for 
HF(v  = 1)  self -deactivation,  when  the  v -♦  R,  R -*  v,  and  R -*  R,T  rate  coeffi- 
cients presented  here  are  used  in  a rotational  nonequilibrium  computer 
program  to  model  the  laser-induced  infrared  experimental  technique  used  to 
measure  this  rate,  will  be  explained. 

0 5 

The  v -*  v and  v -*  R rate  coefficients  were  found  to  have  a T * tem- 
perature dependence  since  the  cross  sections  for  these  processes  were 
practically  independent  of  the  collision  energy.  The  R -*  R rate  coefficients 
have  a complex  temperature  dependence  since  the  cross  section  for  these 
processes  were  found  to  be  dependent  on  collision  energy. 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operation*  of  The  Aerospace  Corporation  it  conducting 
experimental  and  theoretical  investigation*  neceaaary  for  the  evaluation  and 
application  of  • cl  entitle  advance*  to  new  military  concept*  and  ayatem*.  Ver- 
■atility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
peraonnel  in  dealing  with  the  many  problem*  encountered  in  the  nation1*  rapidly 
developing  apace  and  miaaile  ayatema.  Expertise  in  the  latest  scientific  devel- 
opments ia  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophvsics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physic*  Laboratory:  Atmospheric  reactions  and  atmos- 
pherlc  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena.  Including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  iormi  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  mstals. 

Space  Science*  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  auroraa 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  wave*  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphsre,  ionosphere,  and 
magnetosphere;  the  effect*  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  system*. 
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